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Abstract 



Measurements of the S and masses, mass differences, lifetimes and lifetime 
differences are presented. The sample used is much larger than those used 
previously for such measurements. In addition, the S production rates in Z— > bb 
and Z— i> qq events are compared and the position ^* of the maximum of the ^ 
distribution in Z^ qq events is measured. 
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1 Introduction 

This paper presents measurements of the masses and mean hfetimes of S~ and 
and of their mass and hfetime differences, together with a study of S~'^ production in Z° 

hadronic decays. 

Previous measurements of the mass and mean hfetime suffer from low statistics 
compared to E~ measurements, since they came from bubble chamber or hypcron beam 
experiments with a large asymmetry in the production of S~ and S^. The Particle Data 
Group [1] lists only ~ 80 events used for measurement of the mass and 34 for its 
mean lifetime, compared to ~ 2400 events for the S~ mass and ~ 87000 for its mean 
lifetime. The present analysis uses about 2500 S~ and 2300 S^, with small backgrounds. 
The symmetry in the production of particles and antiparticles in Z° decays makes direct 
measurements of E" and mass and lifetime differences with high precision feasible. A 
non-zero value of cither difference would signal violation of CPT invariancc. 

A comparison of the S production rates in Z° bb and qq events is also pre- 

sented, together with a measurement of the position ^* of the maximum of the distribution 
in ^ = — InXp, where Xp is the fractional S momentum. 

2 The DELPHI detector and event selection 

The DELPHI detector is described elsewhere [2,3]. The detectors most important for 
this analysis are the Vertex Detector (VD), the Inner Detector (ID), the Time Projection 
Chamber (TPC), and the Outer Detector (OD). The VD consists of three concentric 
layers of silicon strip detectors, located at radii of 6 cm, 9 cm and 11 cm. The data 
used here were taken in 1992-1995 inclusive, when the polar angles covered for particles 
crossing all three VD layers were 43° < 6 < 137°, where 6 is given with respect to the z 
axis"f. In 1994 and 1995, the first and third layers had double-sided readout and gave both 
R(j) and z coordinates. The TPC is the main tracking device where charged-particle tracks 
are reconstructed in three dimensions for radii between 29 cm and 122 cm. The ID and 
OD arc two drift chambers located at radii between 12 cm and 28 cm and between 198 cm 
and 206 cm respectively, and provide additional points for the track reconstruction. 

A charged particle was accepted in the analysis if its track length was above 30 cm, 
its momentum above 100 MeV/c, and its relative momentum error below 100%. 

An event was classified as hadronic if it had at least 7 charged particles with momentum 
above 200 MeV/c carrying more than 15 GcV reconstructed energy in total and at least 

3 GeV in each hemisphere defined with respect to the z axis. 

The analysis used 3.25 million reconstructed hadronic decays of the Z, consisting of 
0.67 million from the 1992 run, 0.68 million from 1993, 1.29 miUion from 1994, and 0.61 
million from 1995. 

Simulated events were produced using the JETSET parton shower generator [4], and 
then processed with the DELPHI event simulation program delsim [3] which fully sim- 
ulates all detector effects. For each of the years 1992 to 1994, about 1 million fully 
simulated qq events were analyzed in the same way as the real data, and about 0.6 mil- 
lion for 1995. The total number of simulated events used was thus about 3.6 million, 

t Antiparticles are implicitly included unless explicitly stated otherwise. 

•I-In the standard DELPHI coordinate system, the z axis is along the electron direction, the x axis points towards the 
centre of LEP, and the y axis points upwards. The polar angle to the z axis is called $ and the azimuthal angle around the 
2 axis is called </>; the radial coordinate is i? = \/ -\- . 
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comparable to the number of real events. The number of S and S decays generated in 
the simulation was about 89000. 

3 Analysis 

The E" hyperon was studied by a complete reconstruction of the decay chain 
S~ — i> An~ , where A — > pir^. A similar analysis procedure was used previously for ft~ 

reconstruction [5]. 

All pairs of oppositely-charged particles were tried in a search for A candidates. 
For each such pair, a vertex fit was performed by the standard DELPHI y° search 
algorithm^ [3]. A pair was accepted as a A candidate if the x^-probability of the sec- 
ondary vertex fit exceeded 0.1%, the measured fiight distance from the primary vertex of 
the A candidate in the xy plane exceeded twice its error, and the angle between the mo- 
mentum vector sum of the two tracks and the vector joining the primary and secondary 
vertices was less than 0.1 radians (the loss of signal due to this cut has been shown to 
be negligible). The inclusive A reconstruction efficiency was around 19% [3], including 
the 63.9% branching ratio of A — > pir^ [1]. The invariant mass of the A candidate was 
required to be between 1.105 GeV/c^ and 1.125 GeV/c^. 

One by one, the remaining tracks of charged particles that crossed the A trajectory in 
the xy plane were then combined with the A candidate to form a S~ candidate. All S~ 
were assumed to originate from the beam interaction point, which was calculated event 
by event. 

A constrained fit was performed if: 

• the intersection between the A and the charged particle trajectory was more than 8 
mm away from the main vertex in the xy plane; 

• the A and charged particle trajectories were less than 7 mm apart in the z direction 
at the point of crossing in the xy plane; 

• and the charged particle had an impact parameter with respect to the main vertex 
in the xy plane of at least 0.5 mm. 

The fit used was a general least squares fit with kinematical and geometrical con- 
straints apphed to each S~ candidate. The 16 measured variables in the fit were the five 
parameters of the helix parameterization of each of the three charged particle tracks and 
the z coordinate of the beam interaction point (the x and y coordinates were so precisely 
measured that they could be taken as fixed). The two unmeasured variables were the 
decay radii of the S~ and A. The S~ decay point was then determined from this S~ decay 
radius and the Tr~ trajectory while the A decay point was determined from the point on 
the proton trajectory at the A decay radius. The curved E~ track was not measured, but 
calculated in the fit. 

Four constraints required the momenta of the S~ and A at their decay points to be 
in the same direction as the trajectory joining their production and decay positions, two 
required the other tt" to meet the proton at the A decay radius, and the last (seventh) 
constrained the A mass to its nominal value (1115.684±0.006) MeV/c^ For further details 
concerning the fitting procedure, see [6]. 

The pull distributions of the 16 fitted quantities were all approximately normally 
distributed, with mean within ±0.1 and standard deviation 1 within ±0.1, both for 
data and for the simulated events. 

The following cuts were used to select the and S samples: 

^AV^ consists of two oppositely charged particles originating from a neutral particle decaying in flight. 
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Year 


92 


93 


94 


95 


Ms- in data 
Ms+ in data 
Ms± in data 


1321.60 ±0.17 
1321.70 ±0.18 
1321.65 ±0.13 


1321.25 ±0.16 
1321.49 ±0.14 
1321.37 ±0.11 


1321.45 ±0.10 

1321.46 ±0.12 
1321.45 ±0.08 


1321.50 ±0.16 
1321.19 ±0.18 
1321.36 ±0.12 


Mh--Mh+ in data 


-0.10 ±0.25 


-0.23 ±0.22 


-0.02 ±0.15 


0.31 ±0.24 


M3±-1321.3 in MC 


0.14 ±0.07 


-0.02 ±0.07 


0.06 ±0.07 


0.30 ±0.09 


Corrected Ms-: 
Corrected M3+: 
Corrected Ms±: 


1321.46 ±0.18 
1321.56 ±0.19 
1321.51 ±0.15 


1321.27 ±0.17 
1321.51 ±0.16 
1321.39 ±0.12 


1321.39 ±0.12 

1321.40 ±0.14 
1321.39 ±0.10 


1321.20 ±0.19 
1320.89 ±0.20 
1321.06 ±0.15 



Table 1: S and S mass fit results. Values are in MeV/c^. In the simulated ('MC') 
sample, the generated E" mass was 1321.3 MeV/c; the corresponding mass shifts per 
year are used to correct the mass values found in the data. The errors are statistical only. 

• the x^-probability of the fit had to exceed 1%; 

• the S momentum, ps, had to fulfill 1.2 < ^ < 4.2 where ^ = —Inxp and Xp — 
Ps/Pbeam'i this corrcspouds to 0.015 < Xp < 0.3 or 0.7 < ps < 14 GeV/c; 

• the S momentum had to point into the barrel region of the detector (| cos^ |< 0.85); 

• the decay radius of the S in the xy plane had to exceed 2 cm; 

• the decay radius of the S in the xy plane had to be less than the A decay radius. 

Figure 1 shows the right-sign (Att" and Att"*") mass distributions and the S signals 
before and after the cuts were applied. Apart from a difference in mass resolution, the 
agreement between data and simulation was very good. The distributions of the variables 
used in the selection of H candidates are shown in Figure 2 for wrong-sign (Att"^ and Att") 
as well as for right-sign (Avr" and Att^) combinations. 

The fit gave a narrow mass peak from S decays on a small background, as shown 
in Figure 3a; 2478 ± 68 S~ and 2256 ± 63 decays were reconstructed, as shown in 
Figures 3b and 3c. The fitted curves consist of a hnear term for the background, and 
two Gaussian distributions with common mean for the signal. The S mass resolution 
depends on momentum. Therefore the signal is, in principle, the sum of an infinite 
number of Gaussians. But two give a reasonably good fit. The fitted widths of the two 
Gaussians were (2.0 ± 0.1) MeV/c^ and (5.6 ± 0.4) MeV/c^, with a relative fraction of 
1.29±0.18. The corresponding widths from fitting simulated data were (1.8±0.1) MeV/c^ 
and (5.5 ± 0.5) MeV/c^, with a relative fraction of 2.01 ± 0.27. This parameterization of 
signal and background was used in the determination of the S~ and masses. 

The only possible physical background is the decay fl^ — > AK^. The number of fl~ 
reconstructed in our S~analysis is estimated to be at most five, and consequently to have 
no significant infiuence. 

3.1 Measurement of S~ and masses and mass difference 

Table 1 gives the fitted mass and mass difference values for the real data. As already 
described, the signal (sec Figure 3) was represented by two Gaussian distributions with 
common mean and the background by a linear term. 
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Year 


92 




93 




94 




95 




Mxo offset in data 


-0.87 


± 


0.06 


-1.09 


± 


0.05 


-0.75 


± 


0.06 


-0.80 


± 


0.06 


Mko offset in MC 


0.01 


± 


0.05 


0.56 


± 


0.04 


0.38 


± 


0.04 


0.68 


± 


0.04 


Corrected M^o offset 


-0.88 


± 


0.09 


-1.65 


± 


0.08 


-1.13 


± 


0.09 


-1.48 


± 


0.09 


Ma offset in data 


0.14 


± 


0.05 


-0.14 


± 


0.05 


-0.09 


± 


0.06 


-0.07 


± 


0.06 


Ma offset in MC 


0.09 


± 


0.06 


0.01 


± 


0.04 


0.09 


± 


0.04 


0.17 


± 


0.04 


Corrected Ma offset 


0.04 


± 


0.08 


-0.15 


± 


0.07 


-0.18 


± 


0.08 


-0.24 


± 


0.08 


Calculated Ms offset 


-0.13 


± 


0.09 


-0.44 


± 


0.07 


-0.36 


± 


0.06 


-0.47 


± 


0.07 


Resulting M=± 


1321.64 


± 


0.17 


1321.82 


± 


0.14 


1321.75 


± 


0.12 


1321.53 


± 


0.17 



Table 2: Measured offsets from the nominal K° and A masses in MeV/c^, and the corre- 
sponding offsets in H mass, together with the final, resulting S"*" mass values. The errors 
of the corrected K° and A masses include the spread from the simulation smearing, see 
text. 

In order to correct for any bias due to the data processing or to the analysis and 
fit procedure, the mass values obtained from the data were corrected by the difference 
between the values obtained in the same way from the simulated events and the input 
value used in the simulation (1321.3 MeV/c^). As no effect could be identified that might 
affect the S~and S masses differently, the correction was calculated once for each year, 
using the corresponding fully simulated 2=*= sample. Table 1 also shows these corrections, 
and the corrected mass values. The statistical errors of the corrected values contain the 
statistical errors of the simulation. 

The 5=*= mass value averaged over all years was (1321.45 ± 0.05) MeV/c^ with a 
probability for the combination of 33% before correction, and (1321.35 ± 0.06) MeV/c^ 
with a probability of 17% after correction. Thus the average correction amounted to 
(-0.10 ±0.04) MeV/c^ 

3.1.1 Mass scale calibration 

The mass scale was calibrated by determining the A and K° masses in the same way, 
and comparing the resulting values with the known values [1]. The A and K° samples 
used for this purpose were spread over each whole year and their sizes were restricted to 
make it possible to use the same signal and background parameterizations as for the S^. 

The A and K° decays were reconstructed by considering all pairs of oppositely charged 
particles, and the vertex defined by each pair was determined by minimizing the of the 
extrapolated tracks. Consequently, this was a purely geometrical vertex fit, as opposed 
to the mass-constrained fit described above for the S candidates. The measured A and 

mass offsets from their nominal values are shown in Table 2. 

The widths of the A and K° mass distributions arc somewhat larger for data than for 
simulation. A study was made in which the reconstructed variables in simulation, one 
by one, were artificially "smeared" and a corresponding extra measurement error added, 
such that the width of the mass peak in simulation agreed with that in the data. The 
spreads of the shifts obtained by smearing different variables, amounting to 0.05 MeV/c^ 

^'The reconstructed A samples were typically twice as large as the reconstructed S samples while the samples were 
typically 10 times larger 
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for K° and 0.04 MeV/c^ for A, were included in the errors for the corrected offsets quoted 
in Table 2. However, the means of the mass values from the smear ings agreed with the 
"unsmeared" values. 

Offsets of the corrected A and mass values from their known values can arise from 
an error in the correction for dE/dx losses, an error in the assumed magnetic field, or, 
most likely, a combination of the two. 

The S mass offset can be expressed as a function of the K° and A mass offsets, as 

bK dxY^ ( AMk 



where bi are the S, K° and A mass shift coefficients due to magnetic field changes, and 
di are those due to dE/dx correction changes. The values of these coefficients were 
found using Monte Carlo techniques: bs = 0.0805 ± 0.0004, bx = 0.2376 ± 0.0010, 
6a = 0.0438 ± 0.0002, = 0.20 d= 0.02, dx = 0.367 ± 0.017 and dA = 0.162 ± 0.011. 

Inserting the observed K° and A mass shifts and the mass offset coefficients into the 
above equation, taking all errors into account, gave the S mass offsets presented in Table 2. 
The last line of that table gives the final corrected mass values per year. The final 
average corrected mass value is (1321.71 ± 0.06 ± 0.04) MeV/c^ with a probability for 
the combination of 53%. The second error quoted is the unfolded contribution from the 
uncertainty in the calibration offsets. 



3.1.2 Other systematic uncertainties 

The effect of using different parameterizations for the shape of the S mass peak and for 
the background was studied, as well as that of using various fitting techniques (maximum 
likelihood and minimum x^). These variations gave a spread in the final S mass value of 
±0.03 MeV/c2. 

Applying the same "smearing" technique to the simulated S events as was described 
above for the A and K° analysis gave a spread in the final S mass value of ±0.02 MeV/c^. 
Again the average of the values from the smearing study agreed with the unsmeared 
value. 

As a cross-check, the S mass was measured as a function of the momentum of the pion 
from the S~ decay. This is generally the one passing through the most material, and 
it is not affected by the A mass constraint. Thus it is the one most sensitive to dE/dx 
corrections. No systematic effect depending on the pion momentum was observed. The 
S mass was also measured as a function of the polar angle of the S momentum and of 
the observed distance in the xy plane from the beam axis. Again no systematic variation 
was found. 

The total systematic error was thus ±0.05 MeV/c^, as shown in Table 3. 

3.1.3 Results 

The measured average S masses are: 

Ms- = (1321.70 ± 0.08 (stat.) ±0.05 (syst.)) MeV/c^ 
Ms+ = (1321.73 ± 0.08 (stat.) ±0.05 (syst.)) MeV/c^ 
Ms-+s+ = (1321.71 ± 0.06 (stat.) ±0.05 (syst.)) MeV/c^ 

where the systematic errors quoted are common to all three values. 
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Source 


MeV/c^ 


A and K° mass scale 
Fit parameterization 
Simulation smearing 


±0.04 
±0.03 
±0.02 


total 


±0.05 



Table 3: Systematic error contributions to the S mass measurement. 

The systematic errors cancel in the mass difference^, where the small statistical errors 
on the uncorrected values can therefore be fully exploited. The mass difference measured 
in the data is 

Am = Ms- - Ms+ = (-0.03 ± 0.12) MeV/c^ 

which corresponds to a fractional mass difference of 

(Ms- - Ms+)/M„ge = (-2.5 ± 8.7) x 10-^. 

This improves the precision on this CPT violation test quantity by a factor of 3 compared 
to the current PDG value of (11 ± 27) x 10"^ [1]. 

3.2 Measurement of E~ and lifetimes and lifetime difference 

The measurement of the mean lifetimes of the S~ and and their lifetime differences 
uses the S~ and S candidates with a Att invariant mass within ±5 MeV/c^ of the nominal 
mass, where the signal to background ratio is about 6:1. This is the sample for which 
data and simulation were compared in detail in Figure 2. 

The mean lifetimes were estimated using a maximum likelihood fit. The time distribu- 
tion of the combinatorial background was estimated simultaneously in the fit by using the 
wrong-sign combinations. The observed proper time distributions and the fitted functions 
for the wrong-sign and right-sign distributions are shown in Figures 4 and 5 respectively. 
As the mean lifetimes are much shorter for c- and 6-baryons than for a S, all S may safely 
be assumed to originate from the interaction point. 

The proper time was calculated as 

t - dsM^/P^ (1) 

where ds is the fitted fiight distance in the xy plane, Ps is the fitted momentum of the 
S candidate in the xy plane, and Ms is the nominal S mass. 

For right-sign and wrong-sign candidates in the proper time interval 0.04 ns to 2.0 ns, 
the following likelihood function was formed: 

TV^ s s 

^ = nnu) -iiBitj). (2) 

i=i j=i 

The first factor, the F{t) product, represents the right-sign (Avr^, Att^) combinations. 
The second factor, the B{t) product, is an empirical parameterization of the wrong- 
sign (Att^, Att") combinations. The same function B{t) was also used to describe the 

lilt was checked that there was no difTcrencc between masses of A and A, and that the K" mass did not depend on the 
charge of the highest momentum particle in the decay. 
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Year 


92 93 


94 


95 




0.131 ±0.012 0.179 ±0.014 


0.199 ±0.015 


0.167 ±0.016 




0.165 ±0.015 0.146 ±0.013 


0.205 ±0.015 


0.169 ±0.020 


At = Ts- — rs+ 


-0.034 ±0.020 +0.033 ±0.019 


-0.006 ±0.021 


-0.002 ± 0.025 



Table 4: Fit results and statistical errors for S and S lifetime fits. Values are in 
nanoseconds. 



background in the right-sign sample. Thus, by maximizing the joint likelihood function 
C, the background contribution in the right-sign sample was naturally constrained to the 
shape of the wrong-sign distribution. 

The right-sign function F{t) was given by 

m = -^{aoS{t) + B{t)) (3) 

(7^0 + 1 

where S{t) is a normalized probability density function for the observed signal, i.e. it 
is proportional to e(t)e~*/^=, where e(i) is an empirical efficiency parameterization of the 

time-dependent form e'^'^'i+'^^*) determined from the simulation. The relative normalization 
of the signal S{t) and background B(t) in the right-sign sample, ctq, was fixed by the 
observed number of right-sign {Nrs) and wrong-sign (N^jg) events in the fitted time interval 
0.04 ns to 2.0 ns, ctq = ^r^^. 

The background function B(t) was given by 



where 



-1 



and A/i and J^2 are normalization constants for the two F-distributions b(t; (Ji). The value 
/3 = 3 provided a good description of the wrong-sign distribution. The parameters 6i, cti 
and (72 were fitted to the data, together with ts. The fit results for each year are given 
in Table 4. 

The measured S~ and lifetimes are: 

Ts- = (0.165 ± 0.007 (stat.) ±0.012 (syst.)) ns 
rs+ = (0.170 ±0.008 (stat.) ±0.012 (syst.)) ns 
rs-+s+ = (0.167 ±0.006 (stat.) ±0.012 (syst.)) ns 

where the results were achieved by performing the same analysis on the four separate 
years. The hfetimes were taken as the weighted average of the four years. 

In order to minimize the effect of statistical fluctuations, the combined S~ and 
sample was used to evaluate the systematic errors. The following sources of systematic 
errors were considered: 

• the effect on the lifetime fits of the uncertainty in the parameters Ci and C2 in the 
efficiency parameterization was estimated using the simulation to be ±0.005 ns; 
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• the difference between the input and reconstructed hfetimes in the simulation was 
(0.002±0.004) ns; 

• changing the fit range between 0.04 ns and 0.08 ns for the lower boundary and 
between 0.6 ns and 2.0 ns for the upper boundary changed the fitted lifetime by 
±0.004 ns; 

• changing the value of j3 (Equation 5) in the range 1.5 to 4.0 had no significant 
influence on the final results; 

• the of the combination of the four different years was 14.1 for 3 degrees of freedom 
corresponding to a probability of only 0.3%; applying the PDG scaling procedure to 
the combination gives an additional systematic error of ±0.011 ns. 

Note that a S could be reconstructed only if the (anti)proton from the (anti)lambda 
was seen in the TPC Therefore the fact that more antibaryons than baryons interacted in 
the material before the TPC reduced the relative number of reconstructed by about 
10%, but had no significant effect on their lifetime distribution. Thus the systematic 
errors quoted above are common to all three lifetime values. 

The systematic errors cancel in the measurement of the lifetime difference, averaged 
over the years: 

< Ar >=< th- - rs+ >= (-0.002 ± 0.011) ns, 
which gives a fractional lifetime difference of 

(th- - T^+)lTa,erage = "0.01 ± 0.07. 

This quantity, which would indicate violation of CPT invariance if different from zero, 
has a much smaller error than the PDG value of 0.02 ± 0.18 [1]. 

The value of At may also be used together with the world average for the hfetime, 
^PDG _ (0.1639 ± 0.0015) ns, to make a new precise determination of the hfetime 
alone: 

rs+ = t|P^ - At = (0.166 ± 0.011) ns. 

3.3 Measurements of S~ and production 

The parameterization of the signal used in the S mass determination was not used 
to evaluate the efficiencies and production rates, since the broader Gaussian tended to 
become unreasonably wide if left free when fitting substantially smaller data samples. 
Instead, as in the lifetime analysis just described, a fixed interval, this time ±10 MeV/c^ 
around the nominal S~ mass, was used as signal region. The background was estimated 
from the wrong-sign invariant mass distributions. 

The efficiency was determined from simulation and depended on the 5 momentum 
(see Table 5 and Figure 6). The average efficiency was found to be (6.76 ± 0.27 (stat.))% 
for the combined S~ and reconstruction, including cuts and the 63.9% branching 
ratio for A — > p7r~. The error comes from the finite number of simulated events. As 
mentioned earlier, the reconstruction efficiency was ~ 10% lower for than for S^, due 
to differences in the cross-sections for hadronic interactions of particles and antiparticles 
in the detector material. 

All E~ candidates satisfying the standard E~ cuts, and with a Att" invariant mass 
within ±10 MeV/c^ of the nominal E' mass, were considered. The background contribu- 
tion was estimated from the wrong-sign combinations and was subtracted. The measured 
distribution in ^ = — InXp is shown in Figure 6 and Table 5. The < 5~±5^ > production 
rate in the ^ interval 1.4 < ^ < 4.0 was found to be 
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< 5" + 5 >gg= 0.0197 ± 0.0007 (stat.) 

where the statistical error includes the contributions from data and simulation. Extrap- 
olating to the full momentum range using the JETSET prediction gave 

< 5- + 5+ ^ 0.0247 ± 0.0009 (stat.) ± 0.0025 (syst.) 

in hadronic Z decays. This result agrees with the previous DELPHI value of 0.0250 ± 
0.0009 ± 0.0021 [7], obtained using a somewhat different S reconstruction procedure, and 
with the OPAL value of 0.0259 ±0.0004 ±0.0009 [8]. For comparison, the DELPHI tuned 
JETSET 7.3 [9] gives 0.0251 and jetset 7.4 with default parameters gives 0.0273, whereas 
HERWIG 5.9 [10] gives 0.0730. 

The systematic error quoted above has the following two sources. Firstly, according 
to the simulation, 20% of the S~ and were produced outside the range 1.4 < ^ < 4.0. 
An error of 50% of this number was taken as a contribution to the total systematic 
error. Secondly, adding a cut on lifetime, ts > 0.1 ns, and requiring the A candidate 
to be tagged as a 'tight' A (xy flight distance above four standard deviations and 
probability larger than 1%) by the reconstruction program gives a very clean sample. 
The production rate calculated with this sample, extrapolated to the full momentum 
range, is 0.0258 ± 0.0012(stat.), which is the same as that above within errors. The 
half-difference of the rates calculated with the two different sets of cuts was added in 
quadrature to give the total systematic error. The effect of varying the width of the 
signal region was very small. 

From a Gaussian fit in the interval 1.4 < ^ < 4.0, the ^ distribution was found to have 
a maximum at 

= 2.50 ± 0.06 (stat.) ± 0.04 (syst.) 

where the systematic error was evaluated by varying the fit region and the S mass window. 
The JETSET model, with parameters tuned as in [9], gave Qetset = 2.522 ± 0.004(stat.) 
from a similar fit. These values are slightly lower than the OPAL measurement of ^* = 
2.72 ± 0.13 [8]. 

The large statistics of the jetset S~ sample clearly showed that the generated ^ 
distribution was not Gaussian. However, fitting a modified Gaussian form [11] to the 
generated ^ spectrum gave ^* = 2.506 ± 0.004, very close to the result of fitting the 
unmodified Gaussian. Fitting the same modified form to the data, keeping the skewness 
and kiirtosis parameters fixed to the values found in the simulation, gave ^* = 2.51 ± 
0.06(stat). These fits were also to the region 1.4 < ^ < 4.0. 

The above procedure for finding S~ was also applied to Z ^ bb decays. The bb 
events were selected with a lifetime tag algorithm [3,12]. This technique is based on 
the measurement of the impact parameter of each particle relative to the Z° production 
point. Decay products from particles with relatively long lifetimes, like S-hadrons, will 
have large impact parameters. Particles produced in the primary interaction will have 
impact parameters with a spread around zero according to the spatial resolution of the 
detector. From all tracks with a positive impact parameter in an event, the probability for 
the hypothesis that they all came from a single point was calculated. Events in which this 
probability was below 1% were selected as bb events. The joint efficiency to reconstruct a 
S~ decay and simultaneously tag a bb event with this cut was about 3%, with a bb purity 
of 77%. Using these results the production rate of S~ and in bb events was calculated. 
Taking the weighted average of the results from the four years gives the final rate: 

< S- + >^l= 0.0183 ± 0.0016 (stat.) ± 0.0035 (syst.) 



10 



ivioinciixum 

yKjti V / 




X 


p 




t _ 
? — 




- \nxp 


t/incieiicy 


ixeconbxrucxeQ 


At; /bin/event 


9.21-11.24 





202- 


-0 


246 


1 


4- 


-1 


6 


1.7±0.4 


65±17 





00118±0.00042 


7.54-9.21 





165- 


-0 


202 


1 


6- 


-1 


8 


3.7±0.4 


173±21 





00145±0.00026 


6.17-7.54 





135- 


-0 


165 


1 


8- 


-2 





5.0±0.4 


290±28 





00178±0.00024 


5.05-6.17 





111- 


-0 


135 


2 


0- 


-2 


2 


8.2±0.5 


474±33 





00177±0.00016 


4.14-5.05 





091- 


-0 


111 


2 


2- 


-2 


4 


10.5±0.6 


604±35 





00180±0.00014 


3.39-4.14 





074- 


-0 


091 


2 


4- 


-2 


6 


11.9±0.6 


758±37 





00196±0.00014 


2.77-3.39 





061- 


-0 


074 


2 


6- 


-2 


8 


11.7±0.6 


804±38 





00210±0.00014 


Z.Z l—Z. 1 1 


U 


UdU- 


-u 


Uol 


o 
z 


o 

o- 


-3 





1 Q n_i_n a 







00176±0.00011 


1.86-2.27 





041- 


-0 


050 


3 


0- 


-3 


2 


11.9±0.6 


633±31 





00165±0.00011 


1.52-1.86 





033- 


-0 


041 


3 


2- 


3 


4 


8.9±0.5 


433±27 





00151±0.00012 


1.25-1.52 





027- 


-0 


033 


3 


4- 


-3 


6 


6.7±0.5 


221±20 





00101±0.00012 


1.02-1.25 





022- 


-0 


027 


3 


6- 


-3 


8 


4.5±0.4 


133±16 





00092±0.00014 


0.84-1.02 





018- 


-0 


022 


3 


8- 


-4 





2.5±0.4 


64±12 





00081±0.00020 



Table 5: S efficiency and ^ distribution. 

where the systematic error comes from the momentum extrapolation and the sample 
variation of the four years' data. Different cuts on the probability as well as looser S~ 
selections were also tested. The value of < -|- >ijI changed only within ±0.001. 
The DELPHI tuned JETSET 7.3 [9] gives <E~ + >hb =0.0238 and jetset 7.4 with 
default parameters gives 0.0208, whereas HERWIG 5.9 gives < 5~ + 5^ = 0.0523. 

4 Summary 

About 2500 S~ and 2300 decays have been reconstructed from data collected by 
the DELPHI detector in the years 1992 to 1995. 

From this large sample, direct measurements have been made of the S~ and masses 
and their average and difference: 

(1321.70 ±0.08 (stat.) ±0.05 (syst.)) MeV/c^ 
(1321.73 ±0.08 (stat.) ±0.05 (syst.)) MeV/c^ 
(1321.71 ± 0.06 (stat.) ± 0.05 (syst.)) MeV/c 
(-0.03 ±0.12) MeV/c^ 
(-2.5 ±8.7) X 10"^ 

The masses given by the PDG [1] are Ms- = (1321.34 ± 0.14) MeV/c^ Ms+ = (1321.20 
± 0.33) MeV/c^ and (Ms- - Mh+) / Ma^erage = (11 ± 27) xlO'^ Up to now only small 
samples of are referenced by the PDG. 

The S~ lifetime measurement obtained is consistent with the PDG value of (0.1639 ± 
0.0015) ns but has a much larger error. The hfetime difference obtained: 

At = Ts- - Ts+ = (-0.002 ± 0.011) ns 



Ms- = 

Ms+ = 

Ms-+s+ = 

Ms- - Ms+ = 

(Ms- - Ms+)/M„_e = 
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implies 

(ts- - Ts+)/Wape = -0.01 ± 0.07 

and, using the PDG value for r^- , 

rs+ = (0.166 ±0.011) ns, 

which is more precise than our direct measurement, ts+ = (0.170 ± 0.008 ± 0.012) ns. 
The present PDG values [1] for the fractional lifetime difference and for the lifetime 

are (ts- — T'b.+ ) /Taverage = 0.02 ± 0.18 and Ts+ = (0.16 ± 0.03) ns respectively. 

Thus this analysis significantly improves the precision on the fractional mass and 
lifetime differences of S~ and 5^, which test CPT invariance, compared to the present 
PDG values. 

The inclusive production rates for S~ plus S in hadronic Z decays and in Z — > 66 
decays were found to be: 

<5- + 5+>^_ ^ 0.0247 ± 0.0009 (stat.) ±0.0025 (syst.) 
< 5" + >65 = 0.0183 ± 0.0016 (stat.) ± 0.0035 (syst.) 

The JETSET predictions agree with these measurements, whereas the HERWIG predictions 
do not. The maximum of the ^ = — hiXp distribution was found to be at 

C = 2.50 ± 0.06 (stat.) ± 0.04 (syst.) 

in hadronic Z decays. 
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Figure 1: The right-sign (Avr^ and Att^) mass distribution with the S (S~ and added) 
signals in different probabihty bins for data and simulation. The data are represented 
by the points with error bars and the simulation by the histograms, which are normalized 
to the same number of entries: (a) shows the S signal without any other cuts applied 
for events with a fit probability below 1%, (b) shows the S signal without any other 
cuts applied for events with a fit probability above 1%, (c) shows the S signal after 
all cuts given in the text were applied for events with a fit probability above 1%. 
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Figure 2: All the variables used in the S selection for candidates in the mass interval 
M= ± 5 MeV/c^. The histograms are from the simulation and the points with error bars 
are the data. The years 1992 to 1995 have all been added, both for data and simulation. 
The simulation histograms are normalized to the data ones. All variables are shown for 
right-sign (Avr" and Avr"*") and wrong-sign (Avr^ and Avr") combinations after all cuts 
have been made: (a,b) probability, (c,d) ^ = — ln(pH/pbeam), (e,f) cosine of the 
polar angle 9 of the H momentum, (g,h) flight distance of the H in the xy plane, (i,j) 
distance in the xy plane between the A and H decay points. 
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Figure 3: 1992-1995 data: (a) the S~ and S sample, (b) the H~ sample, (c) the S 
sample. The points with error bars show the right-sign (Avr", Att^) combinations. The 
wrong-sign (Att"'", Att") combinations are shown as the shaded histograms. The curves 
show the fits to the S mass distributions described in the text (solid hne). 



16 




Figure 4: The observed time distribution in the wrong-sign sample for 1992-1995 data. 
The two lower curves are the 6-functions described in the text. Their sum, used to describe 
the combinatorial background, is also shown. Only events with times larger than 0.04 ns 
were used in the fit. 




Figure 5: The observed time distribution in the right-sign sample for 1992-1995 data. The 
lowest curve is the estimate of the contribution from combinatorial background events, 
obtained by fitting the wrong-sign combinations. The middle curve is the estimate of the 
contribution of S~ and decays. The upper curve represents the fit to the observed 
time distribution, i. e. the sum of the two lower distributions. 
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Figure 6: (a) Efficiency-corrected distribution of ,^ = — \nXp: tlie points witli error bars 
represent tlie measured ^ distribution, a fit to a Gaussian function is superimposed, and 
the JETSET E, spectrum is sliown as tfie solid histogram, (b) The S~ reconstruction 
efficiency as a function of ^. 



